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ABSTRACT 

We present hydrogen emission line profile models of magnctospheric accretion onto 
Classical T Tauri stars. The models are computed under the Sobolev approximation 
using the three-dimensional Monte Carlo radiative-transfer code TORUS. We have cal- 
culated four illustrative models in which the accretion flows are confined to azimuthal 
curtains - a geometry predicted by magneto-hydrodynamical simulations. Properties 
of the line profile variability of our models are discussed, with reference to dynamic 
spectra and cross-correlation images. We find that some gross characteristics of ob- 
served line profile variability are reproduced by our models, although in general the 
level of variability predicted is larger than that observed. We conclude that this exces- 
sive variability probably excludes dynamical simulations that predict accretion flows 
with low degrees of axisymmetry. 

Key words: accretion, accretion discs - line: profiles - radiative transfer - stars: 
circumstcllar matter - stars: magnetic fields - stars: pre-main-sequence 



1 INTRODUCTION 

Classical T Tauri stars (CTTSs) are understood to be low- 
mass pre-main sequence stars accreting from circumstellar 
discs. Their spectra typically show emission lines of hydro- 
gen and He I, and both permitted and forbidden metal lines 
(see e.g. lBertoutil989h. while red-sh ifted absorption at high 
velocities (e.g lEdwards etafl Il994fl suggests infall to the 
stellar surface. The basis for several models of the accre- 
tion process is a structured stellar magnetic field disrupting 
viscous disc transport near the Keplerian co-rotation radius. 
Gas is loaded onto the field lines and freefalls to produce hot 
impact regions at mid or high latitudes. Photospheric ab- 
sorption line depths are diminished in comparison to those 
of similar main-sequence stars because of the veiling by ex- 
cess continuum from the accretion hotspots. Magnet ospheric 
accretion models are reviewed bv lHartmannl (1998). 

The interaction between a stellar magnetic field and 
a disc can change the angular momentum evolution of 
the star. T Tauri stars usually rotate well below their 
break-up velocity despite accreting substantial amounts of 
material from a circumstellar disc. Details of t he mag- 
netospheres predicted by different models (e.g. iKoenigll 
Il99ll. IShu et al] Il994l) will have different effects on the 
stellar rotation period. Accretion from outside the Ke- 
plerian co-rotation radius can dissipate angular momen- 
tum to the disc, but gas from within the co-rotation ra- 
dius will spin up the s tar. The requirement of a disc for 
this control mechanism JCollier Cameron fc Campbelil9"9l 



lArmitaee fc Clarke! Il99f3) could be used to explain mul- 
tiple populations of PMS stellar rotation rates (see e.j 
IChoi fc Herbstlll996t IStassun et ai]ll999L lHartmanrJ 12002 
Most models assume a dipole field configuration of or- 
der lkG, although the finite radius over which the field 
must couple to the disc will introduce some distortion due 
to differential rotation on either side of the co-rotation 
radius. A more complex geometry with open magnetic 
field lines could launch th e outflows (e.g. IShu et all Il994l : 
iPaatz fc Camenzindlll996l) th at are inferred from obse rva- 
tions of blueshifted absorption. I Johns-Krull et alJ lll999l) de- 
tected a 2kG longitudinal field component from a He I emis- 
sion line in the spectrum of BP Tau, providing some evidence 
for or dered magnetic fields around CTTSs. Furth er observa- 
tions llValenti et al l200aiSvmington et, all2 0051 found that 
the field strength of CTTS varied smoothly over consecutive 
nights and might be modulated by rotation of the system. 

iHartmann et alJ lll994l) presented the first radiative 
transfer models of CTTSs with accretion occuring via an 
axisymmetric dipolar magnetic field structure, and it was 
found that the simulated hydrogen line profiles were suc- 
cessful in reproducing the gross characteris tics of the obser- 
vation s. I mproved models were pr esented bv lMuzerolle et al] 
Jl998l) and lMuzerolle et alJ (2001), incorporating line damp- 
ing and extending the statistical equilibrium calculation of 
the hydrogen atom to 20 levels. These developments led 
to better quantitative agreement with the observed line 
profiles, to the extent that the model has subsequently 
been used to derive mass-accretion rates from Ha? ob- 
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servations (e.g. [Muzerolle et alJ l2003t : iBouvier et alJ Eoolj; 
lLawson et al.ll2004ft . 

Despite these successes, it is clear that the circumstel- 
lar geometry around a typical CTTS is far from a static, 
axisymmetric, dipolar infl ow. For example, monitoring of 
six CTTSs over five nights (ISmith et al.lll999l) revealed per- 
mitted line profile variability attributed to rotation (time- 
scales of days) , variable accretion (several hours) , and flaring 
(about an hour). It was found that the H7 and H<5 line pro- 
file variability was strongly correlated but that there was 
a time delay between variations in the high-excitation lines 
and those fro m lower-energy lines. Hourly variations seen 
in RU Lup by IStempels fc Piskunovl 1)2002;) are attributed 
to varying mass accretion rates and/or changes in the dis- 
tribution of accreting matter. Those authors conclude that 
the observations rule out an axisymmetric accretion flow 
from a disc aligned perpendicular to the stellar rotation 
axis. A spectroscopi c and photometric study of AA Tau 
iBouvier et alJ l2003) revealed a time delay between varia- 
tions in the Ha and H/3 profiles and the veiling continuum. 
It was argued that this lag was the result of a change in 
accretion rate propagating (in freefall) along the field lines, 
with the change in density being observed first in the line 
profiles and then in the continuum when the increased ki- 
netic energy of the flow is liberated at the hot spots. 

Variability is ofte n observed on rotational time-scales. 
Ijohns fc Basril Jl995bT) obtained optical time-series spectra 
of SU Aur and found variations in the red-shifted absorp- 
tion component of H/3 on the same time-scale as the stellar 
rotation period. A blue-shifted absorption feature consistent 
with an outflow was found to be anti-correlated with the red- 
shifted inflow signature, suggesting a magnetosphere that is 
seen to alternatel y favour inflow then outflow at different 
rotational phases. IPetrov et al.l dl996ll also found evidence 
for simultaneous outflow and infall in SU Aur (Balmer lines 
and He I 5876 A). The possibility of an offset dipolar accre- 
tion structure was considered, but the complex kinematics 
of the magnetically channelled gas and the degeneracy of 
the observational diagnostics make quantitative conclusions 
difficult to achieve. 

Long-term m o nitori ng of CTTS TW Hya 
iAlencar fc Batalhal l2002h revealed evidence for both 
inflow and outflow in the Balmer line profiles, with the 
deepest Ha blueward absorption occuring at phases when 
the veiling continuum (and hence mass accretion rate) is 
greatest. It was suggested that the system is viewed at low 
inclination, although the presence of photometric variability 
argued against a truly face-on orientation. A similar viewing 
angle was used to explai n the redwards peaked Balmer line 
pro files of PR Tau |Alencar et al.ll200ll) . 

lOliveira et alJ <|2000T) performed cross-correlation anal- 
ysis of multiple spectral lines from SU Aur to investigate 
non-coincident changes. Different components of line profiles 
were seen to be time-lagged with respect to each other, which 
is inconsistent with the expectation for a simple geometry. 
It was suggested that the accretion flow had an azimuthal 
component and that the observer's line of sight was sampling 
different different physical environments for line formation. 
Structur ed accretion is also supported by the observations of 
VZ Cha iSmith et al .1200 ll) who report features that persist 
for longer than the likely freefall time from the inner edge 
of a truncated disc. The possibility of confinement to two 



rotating streams was considered, with the observed period 
being only half the stellar rotation period. 

If we are to fully understand the magnetospheric ac- 
cretion process, from its role in rotational breaking to the 
use of line profiles as accretion diagnostics, we must obtain a 
better description of the geometry of the circumstellar flows. 
The observational evidence for structured accretion is over- 
whelming, and it is now appropriate to challenge these data 
with theoretical models. In this paper we present the first 
3D radiative-transfer modelling of azimuthally-structured 
accretion. We note that it is not the aim of this work to 
reproduce in detail the variability of permitted line profiles 
in CTTS, but rather to gain insight into the likely charac- 
teristics of line profile variability that would be produced 
by 'cartoon' geometries that have been invoked by other re- 
searchers to interpret their timeseries data. As a first step we 
adopt a canonical, bipolar reference model, and examine the 
synthetic line profiles produced by simple departures from 
axisymmetry in the magnetosphere. 



2 RADIATIVE TRANSFER MODELLING 

We have ext ended the torus radiative-transfer code 
jHarriesI 12000) in order to simulate hydrogen line profiles 
from pre-main-sequence stars with magnetospheric accretion 
flows. The circumstellar environment of accreting T Tauri 
stars encompasses a large range of length scales. Assum- 
ing that the material is strongly influenced by the magneto- 
sphere, we must consider a volume extending to several stel- 
lar radii from the surface. The addition of a disc (or a large- 
scale outflow) will further enlarge the space in which our 
simulated photons may have significant interactions. Close 
to the stellar surface, the magnetic fields may confine gas 
in narrow columns, leading to small impact regions at the 
photosphere. The liberation of energy by shock heating by 
the infalling gas in these areas may be the source of much 
of the system's luminosity. 

The complexity of the circumstellar environment re- 
quires a careful choice of mapping to the simulation's grid. 
Regular subdivision of our model space (e.g. in spheri- 
cal polar coordinates) leads to many of the cells repre- 
senting regions with little effect on the radiation field. 
By adopting the technique of adaptive mesh refinement 
(AMR) we can tailor the subdivision of the simulation 
space to improve the sampling of the structures it contains 
(see e.g. iKurosawa fc Hillierl 1200 j ISteinacker et al.l 120031 : 
lHarries et alJl2004HKurosawa et alJl2004l) . Our whole cube- 
shaped space is divided into eight, smaller cubic subcells; 
and those subcells are themselves repeatedly divided until 
some criterion is satisifed throughout the grid. For the mod- 
els presented in this paper, we specify the maximum mass 
that any subcell can contain. This rule can be applied effi- 
ciently if we can quickly calculate the density at any loca- 
tion within the grid; and it concentrates many small subcells 
where gas is present, while the empty regions are only sam- 
pled by a few large cells. ISvmingtonl (120051) gives further 
details of our AMR implementation. 
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Figure 1. Hydrogen line profiles for an axisymmetric accretion 
geometry with parameters: M = 10 — 7 Mq yr — 1 , accretion stream 
T max = 7000K, R, = 2R , M* = 0.5M Q , viewing inclination 
= 30°. The accretion stream begins in the disc plane between 
r inncr = 2.2R* and r ou tcr = 3. OR*. Solid line: TORUS 3-D com- 
putation with AMR scheme. Squares: MCH computation (repro- 
duced from http://cfa-www.harvard.edu/cfa/youngstars/) 



2.1 Accretion flow model 

In order to investigate the effect of azimuthal structure on 
the line profiles we must adopt a reference model for the 
circumstellar geometry. Here we use the m odel of an accre- 
tion d isc and stellar magnetic field given bv lHartmann et alJ 
lll994l hereafter HHC), in which a dipole magnetic field in- 
tersects the disc and allows gas to free-fall along the field 
lines to the stellar surface. The pressure of the gas picked 
up from the disc is assumed to cause no distortion of the 
magnetic field, and the innermost field line defines the ra- 
dius where the disc is completely truncated. The density 
and velocity at any point in the flow can by computed di- 
rectly from the analytical formulae of HHC, which are based 
on the conservation of mass and energy. We have not at- 
tempted to include rotation of the magnetosphere in our 
model. This simplification may slightly change the central 
peak of our simulated line profiles (dominated by low veloc- 
ity gas near the disc), bu t should not affect the line wings 
(see lMuzerolle et al.l200lL hereafter MCH). Altho ugh not as 
sophi sticated as using a Shu-type magnetosphere iShu et alJ 
I1994T) . or indeed ad opting the density str ucture from an 
MHD simulation (e.g. lRomanova et all20o ll), this is the sim- 
plest, plausible circumstellar structure, and it enables us to 
easily investigate the effect of departures from axisymmetry 
on the emission lines. 

Upon reaching the stellar surface, shock heating of the 
infalling gas will liberate its kinetic energy. We adopt the 
assumption of HHC that the radiation (typically X-rays) 
is immediately reprocessed and all the accretion energy is 
radiated from the impact region as a blackbody. The geom- 
etry and mass-transfer rate of our accretion flow (both of 
which quantities may vary with time) define the pattern of 
hot-spots over the stellar surface. At each temporal phase in 
our simulations we calculate the kinetic energy of the flow 
directly above a grid of points distributed over the stellar 
surface. Ignoring any heating/cooling time-scale, we spec- 
ify the emitted spectrum of each surface element to be the 
sum of the intrinsic stellar spectrum and a blackbody con- 
tribution from the instantaneous local accretion rate. The 
time- independent accretion models presented in this paper 
create hot-spots with temperature ~ 7250K. 

Perhaps the greatest uncertainty in the magnetospheric 



radiative-transfer modelling to date is the temperature 
structure of the accretion flow. The form of the tempera- 
ture run will have a significant impact on the line source 
functions, and hence the line profiles themselves. To our 
knowledge t he only s tudy t o address this aspect explicitly is 
the work of iMartml l|l996l) . He calculated the temperature 
structure of dipolar accretion flows with the same geometry 
and density as HHC. Heating was found to be dominated 
by adiabatic compression as the magnetic field lines con- 
verge along the path of the flow, while the major coolants 
were bremsstrahlung radiation and line emission from Call 
and Mgn. He found that the temperature of the accretion 
flow grew from the disc, and reached peak temperatures of 
~ 6000 K before the temperature fell as the density (and 
hence the cooling) increases towards the stellar surface. For 
the purposes of this paper, in which we are concentrating 
geometrical effects, we have adopted the temperature pro- 
file of HHC (see fig. 6 in t heir paper), which qualitatively 
matches the IMartinl lll996l) results. We defer studying de- 
partures from this temperature law to future work, but note 
that multi-line spectral analysis of observations should pro- 
vide strong constraints on the magnetospheric temperatures. 

The statis tical equilibrium code is based on STATEQ 
jH^mgsJ_|f^2Sli! an implementation of the method of 
iKlein fc Castor! lll97St) . The radiation field is computed by 
integrating the spectrum from each surface element of the 
star, weighted by the solid angle subtended from the centre 
of a grid cell. We typically take advantage of any symme- 
tries present in our models to increase the computational 
efficiency. For example, for axisymmetric accretion flows, we 
can initially calculate the statistical equilibrium for grid cells 
in a plane, and then rotationally map the solution charac- 
teristics to the other cells. The iterative calculations for all 
of the cells throughout the 3D grid will then proceed quickly 
because they are starting close to the solution they will con- 
verge upon. We compute the first 14 levels of the hydrogen 
atom in order to have sufficient high-n levels to avoid the 
inevitable overpopulation of levels near the continuum im- 
pacting on the upper levels of the lines under consideration. 

2.2 Monte Carlo radiative transfer 

The r adiative transfer algorithm is described in lHarriesI 
(2000) with reference to a spherical polar grid. We have up- 
dated the implementation to use recursive routines so that 
it can be applied to our AMR grid, which is stored in a 
tree structure. The emissivity of each grid cell in the flow 
region is used to weight the probability of that cell being the 
source of a Monte Carlo photon. The stellar photosphere is 
separately represented by a grid of small surface elements. 
Photons originating from the star are weighted by the con- 
tinuum luminosity of each element at the wavelength of the 
spectral line. 

The numerical integration from th e sourc e of a photon 
to the observer is described in lHarriesI (|200ff) . Sampling of 
the conditions along each photon's path is done with a re- 
cursive traversal of the tree structure storing the AMR cells. 
The line transfer is performed under the Sobolev approxima- 
tion, and we neglect non-local coupling terms (which would 
render our 3D calculations computationally intractable). 
Both assumptions lead to a negligible effect in the line wings, 
but may have a significant impact at the line core. Addi- 
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Figure 2. The four accretion models (see Table HI : (a) D150 
(b) D90 (c) D30 (d) A30. The observer views the systems at 
inclinations of 30° (a and d) or 60° (b and c). The intensity is 
measured over a narrow wavelength band around the Ha line 
(6550-6575A). The accretion footprints are obvious as bright 
hotspots at high latitudes. The images are logarithmically-scaled 
over three decades from most intense (black) to least intense 
(white). The opaque accretion discs truncate the view of the ac- 
cretion stream(s) below the disc planes. 



tionally, we do not include the effects of line damping. Our 
model profiles are therefore likely to underestimate the true 
extent of the line breadth for Hcv (see MCH for a discussion 
and comparison of results) but will have little effect on the 
higher order lines. 

Alfen-driven turbulence has been invoked to explain 
the broad, s ymmetric wings that are observed in some 
profiles (e.g. | Decamp li||l98ll: lHartmann et al.l Il982l ; iBasril 
ll99fJ: lAl encar & Basri 2000). We have not included turbu- 
lent broadening in the calculations presented here, since 
if the velocity amplitude of the turbulence is comparable 
with the bulk flow velocity then the Sobolev approximation 
breaks down. As we are interested in the imprint of geo- 
metrical structure on the emission line profiles we do not 
regard the lack of turbulent broadening in our models as a 
significant problem, and furthermore we note that its inclu- 
sion requires co-moving statistical equilibrium calculations 
in three-dimensions, a problem that is computational in- 
tractable at present. 

We choose to present our synthetic spectra normalized 
to the local continuum at each time interval. We do not cor- 
rect for veiling of absorption features or suppression of emis- 
sion features by a changing continuum level. While observa- 
tions _of_comjjaxable sp^ectra can compensate for these effects 
(e.g. lEdwards et lull 1994ft . these corrections introduce their 
own uncertaintie s and many authors present spectra wit hout 
such processing dAlencar et alJ^Ooll : ISmith et al.ll200lft . 



2.3 Code testing 

We have performed numerical tests to ensure the spatial 
resolution of our simulations (typically based on ~ 10 6 grid 
cells) is sufficient to represent the accretion flow geometry 
we are using. In Fig.^we compare Pa/3 and Br7 profiles (se- 
lected to minimize discrepancies due to line damping) com- 
puted in 3D using TORUS with the profiles for the same ge- 
ometry calculated using the code described by MCH. Given 
the slightly different physics underlying the two codes the 
profiles show satisfactory agreement in intensity and mor- 
phology. 

The Monte Carlo radiative transfer data sets presented 
later have each been generated from at least 1 x 10 6 
photon packets. The noise correlation between our spec- 
tral wavelength bins, and our adoption of many variance- 
reduction techniques, prohibits simple estimation of the ef- 
fective signal-to-noise ratio, but the excellent reproducibility 
of our results suggests noise is insignificant. 

2.4 Accretion models 

Time-series observations of the line profile variability of 
CTTSs have not yet provided unambiguous clues to the 
structure of the accretion flows. The observer is likely seeing 
a geometry with overall characteristics persisting for longer 
than the rotational time-scale, but with some structure vary- 
ing more rapidly. The accretion rate is also expected to 
show changes on short time-scal es dependent on th e chang- 
ing magnetic field configuration iSmith et alJ 2001). 

In this paper, we aim to identify some of the charac- 
teristics of the magnetically-controlled accretion. If, in one 
simulation, we depart in too many ways from the simple 
steady-state models it will be difficult to identify the con- 
tributions from each change. We have therefore made only 
small adaptations to the HHC magnetospheric accretion ge- 
ometry, breaking the axisymmetry to investigate line pro- 
files as the star rotates. The accretion flows are confined to 
two azimuthal curtains around the star so that their projec- 
tion against the stellar disc and hot-spots chan ges during a 
rotati on cycle. This scenerio was suggested bv lSmith et alJ 
(2001) as an explanation for the periodicity in the spectra 
of VZ Cha. 

This configuration also is similar to that f ound in the 
3D MHD simulations of iRomanova et all ll2003ft . Those au- 
thors model a magnetic field slightly inclined relative to the 
star's rotation axis and the disc plane. Accretion is energeti- 
cally favourable at some azimuth positions, and the flow be- 
comes concentrated in two streams (see fig. 5 of their paper). 
We have simulated accretion curtains which are symmetric 
about the disc plane, and also streams that are diametrically 
opp osite about the star. 

IRomanova et alJ 112003ft found that the accretion was 
axisymmetric in the case where the magnetic and rotation 
axes were the same. For one of our non-axisymmetric mod- 
els, we have therefore implemented an offset between these 
two axes. The HHC formulation for the accretion flow would 
not represent such a misalignment without substantial re- 
working. We have instead used the HHC accretion flow un- 
modified, and achieved the offset effect by only changing the 
rotation axis for our simulated observations, and the disc 
axis. The start of the accretion flow is then not coincident 
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Figure 3. Intensity images of accretion models at selected phases. Each image shows flux integrated over a narrow wavelength range 
(corresponding to ±500kms — 1 ) around the rest wavelength of the indicated spectral line, (a) D90 viewed at the wavelength of Br7 at 
phase (60° inclination). The two streams lie directly in front and behind the star, (b) A view of D90 at phase 0.25 (H7 passband, 
60° inclination). The streams are perpendicular to the line of sight, (c) A30 at Pa/3 at phase 0.5 (30° inclination). The lower stream is 
towards the observer but viewed through the disc's inner hole and so appears truncated. The upper stream is on the far side of the 
star — the accretion hotspot is marginally visible at the limb. The images are logarithmically scaled over four decades, from black (most 
intense) to white (least intense). 

Table 1. Model parameters and properties of the spectral continua. Double accretion curtains are mirrored about the disc plane, while 
the alternate configuration has one stream above the disc, and one below, at opposing azimuth angles. These models are illustrated in 
Fig. |21 The mean levels are quoted for an observer at a distance of 100 pc. The range is the full range of the continuum variability, 
expressed in magnitudes. The final column gives the relative contribution to the continuum flux from the accretion hotspot (f aC c), relative 
to the photospheric flux (f p hot) a t the line wavelength. 



Model 


Accretion curtains 
Type Angular size 


M 

(Moyr- 1 ) 


Magnetic offset 


Line 


Inc. 


Mean continuum 
(ergs - 1 cm _2 A 1 ) 


Range 
(mag) 


face /fphot 


D150 


Double 


150° 


8.3 x 10~ 8 


0° 


H7 


60° 


3.9 x 10~ 13 


0.06 


3.9 


D90 


Double 


90° 


5.0 x 10" 8 


0° 


H7 


30° 


3.7 x 10~ 13 


0.0 


2.3 














60° 


2.7 x 10" 13 


0.15 














Ha 


60° 


3.4 x 10~ 13 


0.07 


0.58 












Pa/3 


60° 


1.7 x 10~ 13 


0.03 


0.16 












Br7 


60° 


5.7 x 10" 14 


0.002 


0.08 


D30 


Double 


30° 


1.7 x 10" 8 


0° 


H7 


60° 


1.5 x 10" 13 


0.14 


0.78 


A30 


Alternate 


30° 


8.3 x 10" 9 


10° 


H7 


30° 


1.3 x 10" 13 


0.3 


0.38 












Ha 


30° 


2.5 x 10~ 13 


0.13 


0.096 












Pa/3 


30° 


1.5 x 10~ 13 


0.04 


0.026 












Br7 


30° 


5.3 x 10~ 14 


0.03 


0.013 



with the surface of the disc, but the results are approxi- 
mately valid for small angular offsets (we choose 10°). 



We adopt the CTTS system parameters of HHC: M, = 
O.8M0, R, = 2R0; and an accretion stream beginning close 
to the disc plane between ri mlcr = 2.2R, and r ou t C r = 3.0R». 
The stellar continuum flux is describe d by a model atmo- 
sphere for T cff = 4000K, log g = 3.5 llKuruczlll993T) . The 
mass accretion rate per unit azimuth angle, per hemisphere 
is the same in each case (integrating to 1O _7 M0 yr -1 if the 
accretion streams were complete), but the limited extent of 
the flow in each model reduces the effective accretion rate. 
Table Q shows the parameters that differ between the mod- 
els. We have simulated line profiles at viewing inclinations 
measured with respect to the stellar rotation axis. 



3 RESULTS 

The results we present are grouped to demonstrate the ef- 
fects of changing individual aspects of the simulation (spec- 
tral line, accretion geometry and viewing inclination). For 
two of the models, we show time-series spectra simulated for 
four spectral lines (Sections 13. II and 13.41 . We also compare 
the spectra for a single line (H7) as the accretion geome- 
try or viewing inclination are altered (Sections 13.21 and 13.31 
respectively). 



3.1 Spectral line comparison 

The D90 model has accretion streams on opposite sides of 
the star, both above and below the disc, i.e. four in total 
(Figs. |5p, 3a, 3b). Each stream covers an azimuth range of 
90°, and the whole system is viewed from an inclination of 
60°. At this angle, the upper accretion hotspots alternately 
rotate out of the line of sight, modulating the oberved con- 
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Figure 4. Spectra from D90 model, viewed at 60° inclination. The bottom plots show the mean spectrum for each hydrogen line. The 
greyscale image shows the quotient of the spectrum at each rotational phase, with respect to the mean spectrum. The variability shown in 
the greyscale has ranges (from top left) of 57, 79, 47 and 35 per cent of the mean spectrum in a given wavelength bin. The plot beside each 
greyscale image shows the continuum level at each phase - the scales are presented in Table HI and the values increase towards the right 
hand side. The RMS deviations of each wavelength bin are plotted in the uppermost box, with the zero point at the bottom of the box. 



timium intensity (Fig. 0] Table 0. This effect is most pro- 
nounced (0.15 mag) at the wavelength of the H/y line because 
the hotspots contribute most significantly to the blue con- 
tinuum. 

The hydrogen line profiles show periodic variations as 
the projected velocity of the inflowing streams changes over 
a rotational period (Figs. 01|5Ji). The mean spectral profiles 
are all peaked at the rest wavelength of the line and are 



asymmetric with most of the emission occuring blueward of 
line centre. Hcv is the strongest of the lines and is in emission 
over its full width; the others all show inverse P Cygni pro- 
files, with H7 having the most significant red-shifted absorp- 
tion feature. The origins of the line profile shapes are found 
in the line source functions (Fig.|SJl. The source functions of 
all four lines lie below their local Planck functions, although 
the lines become more thermalized as the density increases 
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rapidly close to the stellar surface. For H7 and the near in- 
frared transitions, the line source functions are significantly 
smaller than the Planck function for a T ~ 7250K contin- 
uum source: an IPC profile therefore results when the hot 
spots are viewed through the intervening accretion streams. 
The Ha source function is significantly more thermalized, 
and does not show red-shifted absorption. 

Observers often present the temporal var iance spectrum 
(TVS) of their line profil e timeseries (e.g. Ijohns fc Basril 
ll995al:lQliveira et al.ll200ol) . and it is a useful diagnostic for 
quantifying both the level and velocity distribution of line 
profile variability. In Fig. 0]we present the root-mean-square 
(RMS) spectra of our models, which is a comparable diag- 
nostic to the TVS. It is immediately apparent that the Ha 
model shows the strongest variability, and that all lines show 
a multiply-peaked RMS spectrum with the highest variabil- 
ity occurring at line centre. 

The H7 profile varies at three velocity ranges, two of 
which are well correlated with each other (seen in the RMS 
plot as bands centred at 200 and Okms -1 ); and the other 
is anticorrelated with those two (seen at ~ — 130kms _1 ). 
There are two regions between these bands which are almost 
invariant because the integrated line emission seen at some 
velocities remains constant throughout the rotation period. 
Ha shows the same zero- velocity and blueshifted bands, but 
the high- velocity red-shifted band is constant (Ha does not 
show an IPC absorption feature). In the two infrared lines, 
the red-shifted and blue-shifted emission is less variable than 
for H7, but the line centre still shows large changes in in- 
tensity. The maximum of the low- velocity component occurs 
at phases 0.25 and 0.75 when the streams are perpendicu- 
lar to the line of sight, and the projected velocity of the 
inflowing gas is minimized (Fig. At these phases, all 
the line profiles (Figs. 0] are composed of a symmetric 
emission peak around krns" 1 , and (excepting Ha) a small 
red-shifted absorption feature (only a very small volume of 
gas is seen in front of the hot-spot with this configuration). 

Phases and 0.5 (Figs. EKiEt) have the streams aligned 
with the observer, and the projected velocities of the gas are 
maximized. Considering first H7, the stream in front of the 
star is responsible for strong absorption at ~ +220 km s -1 . 
Blueward of line centre is a broad flat-topped emission fea- 
ture extending to ~ — 170kms . The line profile variability 
in the red line wing is strongest for H7, since the hot spot 
emission (which is dominant at short wavelengths) is ab- 
sorbed by the infalling material. The red-wing variability is 
considerably smaller for the other transitions, mainly due to 
the smaller contribution of the hot-spots to the total con- 
tinuum flux. 

3.2 Geometry comparison 

Our second group of results (Fig. [7J are all for the H7 line 
observed from an inclination of 60° , and we show synthetic 
spectra for accretion streams with different azimuthal ex- 
tents. The mean spectrum for the D90 model is shown with 
the other two mean spectra (after rescaling to equal equiv- 
alent widths) for comparison, and the lines have almost the 
same profile. The range of variability in the continuum is di- 
rectly related to the extent of the accretion footprints, with 
the models with the most localized accretion showing the 
most significant continuum modulation. 



(a) (b) 
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Figure 5. H7 spectra at selected phases, (a) D90 model (viewed 
at i = 60°) when curtains perpendicular to line of sight (phase 
0.25; solid line) and along line of sight (phase 0; dashed line), (b) 
A30 model (viewed at i = 60°) with curtains perpendicular to 
line of sight (phase 0.25 or 0.75; solid line), with upper stream on 
far side of star (phase 0.5; dashed line) and with upper stream 
towards observer (phase 0; dot-dashed line). 
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Figure 6. Source functions, S, for spectral lines in the accretion 
columns. Solid line: H7; dot-dashed: Ha; dashed: Pa/3; dotted: 
Br7. The local Planck function B v is shown in the same style. 
For a given spectral line, the upper line shows B v and the lower 
line shows S. 



Although the mean spectrum can not be used to dif- 
ferentiate between different models, the time-series quotient 
spectra could be used as an observational diagnostic. The 
periodic behaviour of the lines is approximately similar in 
the three simulations, but the differing filling factors of the 
accretion flow produce different variations. When the nar- 
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Figure 7. H7 spectra for three models viewed at 60° inclination, see Fig. 4 for a description. From left to right, D150, D90, D30. The 
variability shown in the greyscale has ranges of (from left to right) 39, 57 and 65 per cent of the mean spectrum in a given wavelength bin. 
The plot of the D90 mean spectrum (solid line) also shows the D150 (dashed line) and D30 (dotted line) mean spectra after they were all 
rescaled to have the same equivalent width. 



rowest accretion flow, D30 (Fig. |2J;)), is seen perpendicular 
to the line of sight (phase 0.25 or 0.75) there is very little 
gas projected against the hot-spots, and the line has only 
a small absorption at ~ 120 km s -1 . At the same phases, 
the D150 model (Fig. creates an absorption feature at 
~ 180 km s -1 because the flow extends further in azimuth 
towards the observer and has components with higher pro- 
jected velocity. The emission peak of the line profile shows 
equivalent behaviour for the same reason. 

With only the accretion filling factor differing between 
the three models, the system's rotation will average out the 
contribution from high and low projected gas velocities seen 
at different phases. It is therefore not surprising that the 
mean spectra differ only by a scaling factor. It is clear that 
observations of both line and continuum variability (that 
adequately sample the stellar rotation period) are needed to 
distinguish between the different geometries we have simu- 
lated here. 

3.3 Inclination comparison 

Fig. |S| shows the D90 model viewed at both 30° and 60° 
inclination. On average, the accretion streams have larger 



velocity components parallel to the disc plane than perpen- 
dicular to it, so the line profile extends to higher veloci- 
ties ( — 180-300 kms -1 ) at 60°, compared to 30° inclination 
( — 120 — 250 kms -1 ). At 30° inclination, the continuum vari- 
ation is comparable with the photon noise, but at 60° one 
of the upper accretion hot spots will not be visible at some 
phases and so the variability becomes significant (0.15 mag). 
Additionally, the two lower hotspots are never visible from 
30°, but at 60° one of them will be seen at some phases. 

The overall characteristics of the time-series quotient 
spectra are the same for both simulations, but the differ- 
ences are clearly seen by inspection of the RMS variability 
spectra. At 60°, the line is varying at three distinct wave- 
length ranges, but a more complex pattern is seen at 30°. 
Some of the additional variability will be caused by the ob- 
server seeing (at some phases) through the inner disc hole to 
the high-velocity gas close to a lower accretion hotspot. In 
both simulations, there are lines-of-sight that view a lower 
(beneath the disc plane) accretion stream through an upper 
accretion stream. The projected velocity differences between 
the two streams is generally greater when viewed from 30° 
inclination, so this is also expected to create additional vari- 
ability. 
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Figure 8. H7 spectra of the D90 model viewed at inclinations of 30° (left) and 60° (right). See Fig. [I]for a description. The variability 
shown in the greyscale has ranges of 35 (left) and 57 (right) per cent of the mean spectrum in a given wavelength bin. 



3.4 Offset magnetic axis 

Our A30 model is in spired by the MHD simulations of 
iRomanova et all j2003l) . The accretion streams are inclined 
by 10° with respect to the rotation axis, to represent a sim- 
ilarly inclined dipolar magnetic field (see Section [2.411 . The 
accretion columns subtend 30° in azimuth and, unlike the 
other models, are not symmetric about their midplane. A 
stream exists above the disc plane on one side of the star, 
and below the disc plane on the opposite side (Fig. |5Jl). 
There are therefore only two hotspots on the stellar surface 
instead of the four that are present in the other models. 

The line profile variability through a rotation period re- 
flects the dramatic changes in the observer's view of the gas 
flow (Figs. [5] IHb)- The viewing angle of 30° (with respect to 
the disc normal and stellar rotation axes) keeps the upper 
accretion hotspot visible at all phases, but each gas flow is 
seen alternately towards, then away from the line of sight. 
Additionally, the lower stream suffers from varying occulta- 
tion by the opaque accretion disc (affecting the low velocity 
component) and the stellar disc (blocking the high velocity 
component close to the obscured lower hotspot). 

At phase 0.5 the upper stream is on the far side of the 
star and the lower stream is at maximum visibilty through 



the disc's inner hole (Fig. |^t). All the hydrogen lines are 
composed of two distinct emission features - one redward 
of the rest wavelength and a larger one blueward. Material 
close to the disc with extremely low velocities is only seen in 
the upper stream and is visible throughout the rotation pe- 
riod so the line centre is not particularly variable. The vary- 
ing presentation of the high- velocity gas in the upper stream 
dominates the quotient spectra of Fig. |2l Phase sees the 
upper stream towards the observer, with its maximum pro- 
jection onto the stellar disc. The blue continuum from the ac- 
cretion hotspot is seen through the highest velocity gas and 
so the deepest wavelength bin in the redshifted absorption 
feature occurs at +210 km s -1 in the H7 line. At Br7, the 
unheated stellar surface contributes most of the continuum 
and the large volume of intermediate velocity gas seen over 
the star creates an absoption minimum at only -l-GOkms^ 1 . 
Ha is not seen in absorption at any phase. There is a gen- 
eral trend for an enhancement of the line at two wavelengths, 
symmetrically about the rest wavelength, which moves al- 
ternately to lower and higher velocities through the rotation 
period. 
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Figure 9. Spectra from A30 model, viewed at inclination of 30°. See Fig.|4]for a description. The variability shown in the greyscale has 
ranges (from top left) of 62, 112, 23 and 12 per cent of the mean spectrum in a given wavelength bin. 



3.5 Correlation of variability 

Observational studies of line profile variability in CTTS fre- 
quently use cross-correlation images to interpret the often 
complex phenomena that occur across profiles and between 
profiles fro m differ e nt tra nsitio ns of species; a t echniq ue pi- 
oneered by Ijohnsl (Il994l) and Ijohns fc Basril jl995af) . We 
have computed auto-correlation (AC) images (see Figs. 1101 
and llll for each line of our mod el time-series, following the 
method of lOliveira et al] i2000l) . The grey border of each 
image corresponds to the continuum, with white regions de- 



scribing positive correlation and black features indicating 
negative correlation. 

The H7 AC image for the D90 model at i = 60° shows a 
region of positive correlation around the rest velocity of the 
line and at high redwards (200 — 250 km s" 1 ) and most blue- 
wards (< 50kms _1 ) velocities. There is a 'cross' of low cor- 
relation centred on (+130, +130 kms^ 1 ) and strong anticor- 
relation between the red side of the profile (0 — 250 kms -1 ) 
and the blue side of the line (< —70 kms -1 ). This anticor- 
relation is intrinsically linked to the diametric nature of the 
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Figure 10. Images representing the auto-correlation function for the H7 line. The accretion model and the viewing inclination are shown 
above each image. The brightest shades represent the strongest correlation; the darkest shades arc the strongest anti-correlations. 
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Figure 11. Images representing the auto-correlation function for some spectral lines. Top row: D90 model at 60° inclination; bottom row: 
A30 at 30° inclination. Within each row, the columns are, from left to right, Ho, Pa/3 and Br7. The brightest shades represent the strongest 
correlation; the darkest shades are the strongest anti-correlations. 
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accretion curtains: when one curtain is in front of the star 
(leading to enhanced absorption), the other is behind the 
star (leading to enhanced emission at blue velocities). The 
'cross' is related to the D90 RMS spectra in Fig. [7| which 
shows a minimum at ~ +130 kms -1 . The discontinuity be- 
tween positive and negative correlation at ~ —70 kms -1 
corresponds to another minimum in the RMS spectra, and 
it can clearly be seen in the grayscales that the intensity 
varies in antiphase on either side of this velocity. 

The D150 and D30 models at 60° inclination have sim- 
ilar features, and so do the D90 models in the near infrared 
Pa/3 and Br7 lines. The D90 Ha line has an anti-correlation 
region between the rest velocity and velocities greater than 
+ 100 tons" ; this additional feature arises because Ha is the 
only line in the D90 i = 60° model that does not show an 
IPC profile. The same effect is seen in the Ha line of the 
A30 model, whose other lines are positively auto-correlated 
over most of their wavelength ranges. 

We have also computed cross-correlation (CC) images 
(with zero lag) of the Ha profiles with the H7 profiles for 
each of our models (Fig. 1121 . For the D90 model the cross- 
correlation image demonstrates that the variability in the 
red sides of the Ha and H7 profiles is anticorrelated. This 
anticorrelation is directly visible in the grayscale images of 
Fig. 0] with bright features (excess emission) on the red side 
of Ha coinciding with increased absorption in H7. A simi- 
lar pattern is seen for the D150 and D30 CC images. The 
A30 CC image shows that the red wing of the Ha profile 
is anticorrelated with practically the whole H7 line. This is 
a consequence of the red wing varying for only fraction of 
the rotation period - the same range of phases at which the 
upper hotspot is visible and suppressing the emission line. 



4 DISCUSSION 

The illustrative models presented in this paper represent a 
first step in quantifying the emission line profile variability 
observed in CTTS. We did not adjust the simulation pa- 
rameters to achieve agreement with any existing data, and 
our choice of simple configurations has produced line profiles 
with high levels of variation. It is therefore most useful to 
consider the general characteristics seen in our time-series 
spectra. 

Our synthetic spectra typically predict levels of line pro- 
file variability that are significantly greater than those ob- 
served. Only the D150 model, which has wide accretion cur- 
tai ns, shows vari a bility that is comparable to published data 
e.g. ISmith et al.| i200ll) . The geometry with the least sym- 
metry (A30) leads to line profiles that change substantially 
in shape and intensity from one rotation phase to another. 
This geometry was selected to repre sent the structures ind i- 
cated by recent MHD calculations IIRomanova et alJ l2003). 
but it appears that the line profiles predicted by such a 
magnetosphere are sufficiently at odds with spectroscopic 
observations that this form of MHD model may be rejected. 

MHD simulations must balance the interactions of the 
accretion disk with a magnetic field that includes large scale 
structure (e.g. a dipolar component) with local irregular- 
ities. Both accretion fl ows onto the star and magnetically 
launched outflows (e.g Mill er fc Stond fl997T) must be de- 
rived self-consistently to explain observed spectral features. 



CTTSs are probably more fully enveloped in accreting gas 
than our A30 model proposes, but equally, the broad sym- 
metric curtains of our D150 and D90 models are too regular 
to represent the circumstellar environment. Better agree- 
ment might have been reached by keeping the high density 
streams, but elsewhere setting a lower mass accretion rate 
that would increase the azimuthal filling factor of the mag- 
netosphere. 

The RMS spectra of our simulations often show several 
distinct minima, corresponding to velocities where the emis- 
sion from line-of-sight iso- velocity surfaces remains constant. 
This pattern of variability is not often seen in the observa- 
tions, altho ugh Ha, H/3, and Na l doublet variance profiles 
of SU Aur llOliveira et alJl2000|) do show multiple regions 
of mininum variability. More typicall y the line variabi lity 
is m ost pronounced eithe r on the red iSmith et al]l200"ll) or 
blue iAlencar et all200lfl sides of the profile. Occassionally a 
bimodal distribution of variability is seen across the emission 
lines, with equal variability in th e red- and blue-line-wing s 
and a minimum at the line centre iAlencar fc Batalhal2002li . 
Given that our models neglect any outflow component, our 
RMS spectra typically indicate strongest variability on the 
red side of the model profiles. However the A30 model, with 
its narrow curtains, shows that blue-side line variability can 
be dominant (particularly at Ha) even with an infall geom- 
etry. 

Cross- and auto-correlation images are widely used 
as a diagnostic of li n e profile variability i n CTTS (e.g. 
[ Johns fc Basril Il995al: lOliyeira et alJ l2000l lAlencar et all 
boo it lAlencar fc Batalhal l^Oojr A square, cross-shaped 
auto-correlation f unction is often seen in the lower-order 
Balmer lines (e. g. iJohn s fc Basrilll995al : lAlencar fc Batalhal 
120021 lOliveira et al]l2000h " and this has been interpreted as 
resulting f rom viewing the the mag netosphere at low in- 
clinations (lAlencar fc Batalhal 1200211 . We find such a pat- 
tern in our AC images for essentially all our models (see 
Figs. 1101 II II , indicating that the AC images may be of lim- 
ited diagnostic potential. The cross-correlation images at 
zero lag (see Fig. I12H show a variety of structures, but we 
find no obvious similarities with the published CC images. 

Interestingly, while the variability of the emission line 
profiles in our simulations is excessive, the ranges of con- 
tinuum flux (Table are comparable with, or even lower 
than, results from photometric surveys, where V band mag- 
nitu des are seen to chang e periodically by up to ~ 1 mag 
(e.g. iBouvier et al.l I1993T) . In the models presented here, 
any smoothing of the azimuthal magnetospheric structure 
(broadening of the curtains for example) in order to re- 
duce line profile variablity, will lead to an equal smooth- 
ing of the surface hot spots (and a corresponding reduction 
in continuum variability). The simultaneous satisfaction of 
the dual spectroscopic and photometric observational con- 
straints represents a challenging problem. 



5 CONCLUSIONS AND FUTURE WORK 

We have presented the first 3D radiative transfer calcula- 
tions of non-axisymmetric accretion onto classical T Tauri 
stars. Hydrogen line profiles throughout the stellar rotation 
period are shown for three systems that have accretion con- 
fined to curtains filling part of a dipolar magnetosphere. 



Emission-line profile modelling of structured T Tauri magneto spheres 13 



D150 i = 60° 



D90 i = 60° 



D30 i = 60° 



A30 i = 30° 






Figure 12. Cross-correlation maps for the He? and H7 spectral lines in the four accretion models. The model and the viewing inclination 
are shown above each image. The brightest shades represent the strongest correlation; the darkest shades are the strongest anti-correlations. 



Three of the models had curtains that were symmetric about 
the disc plane, and the fourth accreted from above the disc 
on one side of the star, and from below on the other side. 
The latter configuration was suggested by MHD simulations 
of stars with an inclined dipolar field, but we find that our 
model creates line profiles that are excessively variable when 
compared to observations. We suggest that CTTS magne- 
tospheres may have high filling factors of accreting gas, but 
with local enhancements that may be in the form of streams 
confined in azimuth. 

The simple models presented here are, of course, un- 
able to reproduce the wide range of variable phenomenon 
observed. Nonetheless, some gross characteristics of the line 
profile variability are reflected in the models, including the 
form of the variability's velocity distribution and the auto- 
correlation images. In the future we will investigate the line 
and continuum diagnostics of non-sta tic accretion models . 
Time-series spectroscopy of AA Tau bv lBouvier et alJ J2003l> 
already hints at changes in the global accretion rate on a 
time-scale of hours. Using a simple freefall model we should 
be able to examine how a change in the mass-accretion rate 
affects the line profiles and continuum flux as it propagates 
down the field lines. Our self-consistent model for the repro- 
cessed accretion power should enable us to determine the 
lag between the line and continuum variability. 

Furthermore, the observation of anticorrelation betwen 
the blue and red side variability of Ha in SU Aur has been 
interpreted using a dipole offset m odel with simultaneou s 
dipolar accretion and an outflow ll Johns fe Basrilll995bl) . 
The inclusion of a wind component in our geometry will 
enable us to investigate such scenarios and make quantative 
predictions regarding the line variability. 

Finally, the TORUS code includes computation of the 
Stokes vectors for the simulated photon packets, and can 
measure linear polarization by scattering events. Addition of 
a magnetic field configuraton to our simulation grid would 
allow circular spectropolarimetry to be simulated. This tech- 
nique has already been used observationally to show CTTS 
emissi on lines originate in str u ctures with a net magn etic 
field Jjohns-Krull et alJ Il999t ISvmington et alJ 120051) . so 
self-consistent modelling of line intensity and polarization 
would allow more stringent tests of magnetospheric accre- 
tion models. 
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